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ON THE SIZE OF THE NON-THERMAL COMPONENT IN THE RADIO 

EMISSION FROM CYG OB2 #5 

Luis F. Rodriguez/ Yolanda Gomez, ^ Laurent Loinard,^ and Amy J. Mioduszewski^ 

RESUMEN 

^^ , Cyg 0B2 #5 es un sistema binario de contacto con emision variable de radio. 

Cn ' Esta emision tiene un estado bajo de flujo dominado por el viento estelar ionizado y 

^-^1 un estado alto de flujo, donde una componente no-termica adicional aparece. Ahora 

^ C^ ' se sabe que las variaciones tienen un periodo de 6.7±0.2 afios. La componente no- 

termica ha sido atribuida a varios agentes: una envolvente en expansion eyectada 
periodicamente por la binaria, emision de una region de choque de vientos, o una 
estrella con emision no-termica en orbita excentrica alrededor de la binaria. La 
determinacion del tamano angular de la componente no-termica es crucial para 
, discriminar entre estas posibilidades. Presentamos un analisis de datos de archivo 

■^ ' del VLA hechos a 8.46 GHz en 1994 (en estado bajo) y 1996 (en estado alto), que 

rn I nos permiten sustraer el efecto de la componente termica persistente y estimar un 

tamaiio angular < 0"02 para la componente no-termica. Este tamafio compacto 
favorece la explicacion en terminos de una estrella con emision no-termica o de una 
region de choque de vientos. 
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tJ ; ABSTRACT 

d , Cyg 0B2 #5 is a contact binary system with variable radio continuum emis- 

sion. This emission has a low-flux state where it is dominated by thermal emission 
from the ionized stellar wind and a high-flux state where an additional non-thermal 
^ ' component appears. The variations are now known to have a period of 6.7±0.2 yr. 

CN ! The non-thermal component has been attributed to different agents: an expand- 

~~J^ ' ing envelope ejected periodically from the binary, emission from a wind-collision 

region, or a star with non-thermal emission in an eccentric orbit around the binary. 
The determination of the angular size of the non-thermal component is crucial to 
discriminate between these alternatives. We present the analysis of VLA archive 
observations made at 8.46 GHz in 1994 (low state) and 1996 (high state), that al- 
low us to subtract the effect of the persistent thermal emission and to estimate an 
angular size of < 0''02 for the non-thermal component. This compact size favors 
the explanation in terms of a star with non-thermal emission or of a wind-collision 
^ ' region. 

C^ ■ Key Words: STARS: INDIVIDUAL (CYG 0B2 #5) — RADIO CONTINUUM: 

STARS 

1. INTRODUCTION 

Cyg 0B2 #5 (V729 Cyg, BD +40 4220) is an eclipsing, contact binary system consisting of two 0-type 
supergiants with a 6.6-day period (Hall 1974; Leung & Schneider 1978; Rauw et al. 1999; Linder et al. 2009). 
As several other luminous 0-star systems in the Cyg 0B2 association, this source was found to show evidence 
of variable radio emission (Persi et al. 1985, 1990; Bieging et al. 1989). The radio emission appears to have 
two states: a low-flux state of ^2 mJy at 4.8 GHz where the spectral index is consistent with thermal emission 
from an ionized stellar wind, and a high-flux state of ~8 mJy at 4.8 GHz where the spectral index is flatter 
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than in the low state. The variations were proposed to have a 7-year period (Miralles et al. 1994) and have 
been attributed to variable non-thermal emission from an expanding envelope arising periodically in the binary 
(Persi et al. 1990; Bieging et al. 1989; Miralles et al. 1994). 

In addition, Abbott et al. (1981) and Miralles et al. (1994) reported on the existence of a radio "companion" 
0.8 arcsec to the NE of the main radio source (which is associated with the eclipsing binary). Observations by 
Contreras et al. (1997) revealed that this radio source has an elongated shape and lies in-between the short- 
period binary and a third star, which was first reported by Herbig (1967). Contreras et al. (1997) suggested 
that the proposed NE radio companion actually corresponds to the wind interaction zone between the binary 
system and the tertiary component. Recently, Kennedy et al. (2010) reanalysed all VLA observations of Cyg 
0B2 #5 and showed that the primary radio source, associated with the eclipsing binary, varies with a period of 
6.7±0.2 yr while the flux from the secondary NE source remains constant in time. These authors proposed that 
the variations in the main radio component can be represented by a simple model in which a fourth star (with 
constant non-thermal emission) moves around the eclipsing binary in an eccentric orbit and the varying radio 
emission results from the variable free-free opacity of the wind between the star and the observer (Kennedy 
et al. 2010). These authors estimate a major axis of 14 AU for the star with non-thermal emission. In this 
scenario, we have a quadruple star system: the contact binary, the star with non-thermal emission in a 6.7 yr 
orbit around the contact binary, and the NE component. 

To advance our understanding of the nature of the time- variable non-thermal emission, a determination of 
the angular size of the source producing it is needed. Using VLA data taken at 8.46 GHz in 1991 October 3 
(when Cyg 0B2 #5 was in the high radio state), Miralles et al. (1994) estimated an angular size of 0''046±0''006 
for the whole (non-thermal plus thermal components) emission. Using MERLIN data taken at 4.8 GHz in 1996 
November 14 (when Cyg 0B2 #5 was again in the high radio state), Kennedy et al. (2010) estimated an 
angular size of ~ 0''077 for the whole (non-thermal plus thermal components) emission associated with the 
primary radio component. Since the thermal contribution comes from the ionized wind that is known to be 
extended, these size estimates can be considered as upper limits to the size of the non-thermal emission. In this 
paper we present the analysis of archive 8.46 GHz continuum observations made with the Very Large Array 
(VLA) of the NRAcHin the A configuration toward Cyg 0B2 #5 in two epochs. During 1994 April 9 (1994.27), 
the source was in the low radio state and only the thermal emission from the ionized stellar wind was present. 
These observations were used to estimate the characteristics of the persistent thermal component. During 
1996 December 28 (1996.99), the source was in the high radio state, with both the thermal and non-thermal 
components present. Using the information obtained from the 1994 data, it is possible to subtract in the (u,v) 
plane the contribution from the thermal component and obtain a more stringent estimate of the angular size 
of the non-thermal emission. 

2. DATA REDUCTION 

The archive data were edited and calibrated using the software package Astronomical Image Processing 
System (AIPS) of NRAO. The parameters of the observations are given in Table 1. The absolute amplitude 
calbrator for both epochs was 1331+305, with an adopted flux density of 5.21 Jy. The (u,v) data were self- 
calibrated in amplitude and phase. An image was made with the 1994.27 (u,v) data and the clean components 
of this image were subtracted from the 1996.99 (u,v) data using the AIPS task UVSUB. The resulting (u,v) 
data, that we refer to as 1996.99-1994.27, was analyzed both in the image and in the (u,v) planes. In the 
image plane a compact, basically unresolved source is observed. A Gaussian ellipsoid fit to this source in 
the image plane using the AIPS task JMFIT indicates maximum possible deconvolved angular dimensions of 
0''035 X 0"029 (major axis x minor axis) for the source. However, White & Becker (1982) have shown that it is 
better to study the angular dimensions of marginally resolved sources directly in the (u,v) plane. We then fitted 
directly the (u,v) data in two dimensions using the AIPS task UVFIT and a Gaussian ellipsoid model. The fit 
gave 0''023 ± 0''013 x 0''023 ± 0''019 (major axis x minor axis) for the source. This result indicates that the 
source is very compact and that, within the error, it shows no strong departures from a circular morphology. 
We then adopt as an estimate of the Gaussian angular size (full width at half maximum) of the source the 
geometric mean of the major and minor axes, that gives 9g = 0''023 ± 0"012. 
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TABLE 1 
ARCHIVE DATA USED 



Epoch 


Project 


Frequency 

(GHz) 


On-source 
Time (min) 


Phase 
Calibrator 


Bootstrapped 
Flux Density (Jy) 


1994 Apr 09 (1994.27) 
1996 Dec 28 (1996.99) 


AR277 
AR277 


8.46 
8.46 


37 

45 


2007+404 
2007+404 


2.84+0.01 
2.89+0.01 



'^Major axis x minor axis; position angle, for a (u,v) weighting of ROBUST — (Briggs 1995). 



The angular size determined for the non-thermal emission depends on the model assumed to fit the (u,v) 
data. A model of a circular disk with constant brightness gives an angular diameter of 9d = 0''019 ± 0''010. 
Given that the errors are comparable to the measurements, we will take this determination as an upper limit, 
9(") < 0''02, to the angular size of the non-thermal emission. 



3. INTERPRETATION 

The small angular size found for the non-thermal component, 9{") < 0''02, can be used to discuss the 
nature of the non-thermal component. An expanding envelope ejected periodically from the binary (Persi et 
al. 1990; Bieging et al. 1989; Miralles et al. 1994) would have to be bigger than the radio photosphere to be 
outside the optically thick region around the contact binary system. We can use the results for the 1994 data 
discussed in the Appendix to estimate that the radio photosphere at 8.46 GHz has an angular size of d{") ~ 72 
milli-arcseconds. We then conclude that an envelope cannot explain the observations since the non-thermal 
emission would be significantly larger than observed. Recently, Pittard (2010) has proposed that the time 
variable emission observed in 0+0 binaries can come from the thermal emission from the wind-collision region 
between the stars. In the case of Cyg 0B2 #5 the wind-collision region between the members of the contact 
binary is probably a strong emitter, but lies deep inside the radio photosphere of the system, although it is 
known that the radio photospheres can be clumpy or asymmetric, allowing the detection of sources embedded 
in it (e.g. WR140; White & Becker 1991). The wind-collision region could also be located between the contact 
binary and a third star. Kennedy et al. (2010) discussed the variations arising from non-thermal emission in 
a wind-collision region between the contact binary and a third star. They demonstrate that thermal emission 
variations that may arise in such a wind-collision region will not have sufficient amplitude to account for the 
variations, and conclude that the thermal-emission Pittard models would predict smaller flux changes than 
observed. The Kennedy et al. (2010) models also account well for the variable radio emission. We then favor 
the possibility of a third star with non-thermal emission in orbit around the contact binary, as proposed by 
Kennedy et al. (2010), or a wind-collision region. 

Even when the non-thermal emission from young, low-mass stars is typically associated with gyrosynchrotron 
emission from Lorentz 7 factors of a few (e.g. Ray et al. 1997), while the emission from massive stars is believed 
to be synchrotron from relativistic electrons with very large 7 factors (e.g. Pittard & Dougherty, 2006), we 
cannot rule out the possibility that the non-thermal emission from Cyg 0B2 ^5 is being produced by a young, 
low-mass star (i.e. a Tauri star) that formed coevally with the rest of the system. The non-thermal emission 
from young, low-mass stars is usually very compact. For the T Tau Sb star, with a mass of 0.6 Mq, Loinard 
et al. (2005; 2007) estimate an angular size of 0.5 mas, that at a distance of 147.6 pc corresponds to a radius 
of ~8 Rq. For the more massive SI star (6 Mq) in Ophiuchus, Loinard et al. (2008) estimate an angular size 
of 0.95 mas, that at a distance of 116.9 pc corresponds to a radius of ^^12 Rq. At the larger distance of 925 pc 
(Linder et al. 2009), we expect a non-thermal star associated with Cyg 0B2 ^5 to have an angular size well 
below the milli-arcsecond. 

Cyg 0B2 ^5 will be again in high state during 2010 and 2011 and a very long baseline intcrferometry study 
is needed to determine accurately the angular size and morphology of the compact non-thermal emission. 
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Fig. 1. (Left) Real (filled squares) and imaginary (empty squares) components of the emission at 8.46 GHz for 1994 
April 9. The dashed line is the least squares fit to the real component. (Right) Same as in the left panel for 1996 
December 28. The imaginary component is consistent with zero, indicating that the source is symmetric about the 
phase center (the origin of the visibility plane) and has no significant structure on these spatial scales. 

4. CONCLUSIONS 

We present the analysis of VLA archive data of the Cyg 0B2 jj^b system taken during the high (1996) 
and low (1994) radio states. This analysis allows the subtraction, in the (u,v) plane of the persistent thermal 
component and a better estimate of the angular dimensions of the non-thermal component. Wc obtain an 
upper limit, 0(") < 0''02 for the angular size of the non-thermal emission. Cyg 0B2 ^5 will be again in high 
state during 2010 and 2011 and a very long baseline interferometry study is needed to determine the angular 
size and morphology of the compact non-thermal emission. 

We are thankful for the support of DGAPA, UNAM, and of CONACyT (Mexico). This research has made 
use of the SIMBAD database, operated at CDS, Strasbourg, France. 

APPENDIX 

A. ONE-DIMENSIONAL ANALYSIS OF THE VISIBILITY FUNCTION 



If the sources have approximate circular symmetry, the analysis of the (u,v) data can be made in a one- 
dimension approximation. To test if this was possible with the data discussed here, the NE component was 
subtracted in two dimensions in the (u,v) plane of both data sets using Gaussian ellipsoid models to leave only 
the emission from the main component, the one directly associated with the contact binary. This subtraction 
was performed using the AIPS task UVMOD. We then fitted the (u,v) data in two dimensions using the AIPS 
task UVFIT and a Gaussian ellipsoid model. The fits gave 0('079 ± 0"007 x 0"066 ± 0"008 (major axis x minor 
axis) for the 1994.27 data and 0"045 ± 0^006 x 0^036 ± 0^005 (major axis x minor axis) for the 1996.99 data. 
These results indicate that, within error, the emission at both epochs does not show strong departures from 
a circular assumption. We also note that at a distance of 925 pc (Linder et al. 2009), the semimajor axis of 
the non-thermal star (7 AU) will subtend an angle of only 0.008" , below the upper limit of 0.020" found for 
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the angular size of the non-thermal component. It is quite possible that this non-thermal star, if present, will 
introduce asymmetries but future VLBI observations will be needed to test this. 

In Figure 1 we present the one-dimensional real and imaginary components of the amplitude as a function 
of baseline. These values were obtained using the AIPS task UVPLT that averages in circular rings around the 
origin of the (u,v) plane. Before applying this task, the source emission was centered in this origin. 

For a marginally resolved wind source, the real part of the visibility, V(b), can be fitted in one-dimension 
with a linear function (Escalante et al. 1989): 

V{b) = 5o - Ab, (1) 

where 5*0 is the flux density at zero spacing (the total flux density), given in niJy, A is the slope of the linear 
fit, and b is the baseline separation in wavelengths. For 1994 only the thermal (T) emission was present and 
we obtain S'o(r) = 3.79 ± 0.09 mJy and A(T) = (1.64 ± 0.18) x 10"*^ mJy wavelengths-^ For 1996 both the 
thermal and non-thermal (T+N) emissions were present and we obtain So(T + N) = 9.46 ± 0.05 mJy and 
A{T + N) = (2.07 ± 0.14) x 10"^ mJy wavelengths"^. The fact that both epochs give similar slopes suggests 
that, as expected, most of the extended emission comes from the thermal component. Subtracting the 1994 fit 
to the 1996 fit we obtain the non-thermal (N) contribution alone, that is given by So{N) = 5.67 ± 0.10 mJy 
and A{N) = (0.43 ± 0.23) x 10^^ mJy wavelengths^^. The characteristic angular size of the emitting region is 
given by (Escalante et al. 1989; Miralles et al. 1994): 

d{") ~ 1.44 X 10^ ^, (2) 

Do 

which gives 9{") ~ 0''011 ± 0''006 for the non-thermal component. It should be noted that our fitting of the 
non-thermal component to a thermal wind model may not be consistent since the structure of the visibilities for 
an optically thin non-thermal source and a partially optically thick stellar wind are not necessarily the same. 
In Table 2 we summarize the parameters obtained for the thermal and thermal plus non-thermal states as well 
as the subtraction, that traces the non-thermal emission alone. 

TABLE 2 

PARAMETERS FROM FITS 





5o 


A 


e 


Emission" 


(mJy) 


(10-^ mJy A-^) 


(mas'") 


T 


3.79±0.09 


1.64±0.18 


62±7 


N+T 


9.46±0.05 


2.07±0.14 


32±2 


N 


5.67±0.10 


0.43±0.23 


11±6 



"T = Thermal, N = Non-thermal, 
'^mas = milli-arcseconds. 
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